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Abstract 

The total cross section of the reaction pp — > ppK'^K~ has been determined 
for proton-proton reactions with Pbeam = 3.67 GeV/c. This represents the first 
cross section measurement of the pp -^ ppK^K^ channel near threshold, and is 
equivalent to the inclusive pp — > ppK~X cross section at this beam momentum. 
The cross section determined at this beam momentum is about a factor 20 lower 
than that for inclusive pp — > ppK~^X meson production at the same CM energy 
above the corresponding threshold. This large difference in the K^ and K~ meson 
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inclusive production cross sections in proton-proton reactions is in strong contrast 
to cross sections measured in sub-threshold heavy ion collisions, which are similar 
in magnitude at the same energy per nucleon below the respective thresholds. 
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1 INTRODUCTION 



Currently there is much interest to determine the total K~ meson production 
cross section in nucleon-nucleon reactions near threshold. This cross section 
is of particular importance in heavy-ion physics. In heavy ion collisions it 
has been observed that the inclusive production cross sections for K^ and 
K^ mesons are nearly equal in sub-threshold reactions when measured at 
the same energy per nucleon below the production thresholds for the reac- 
tions pp — * K~X and pp -^ K^X, respectively [1]. This result is surprising 
for several reasons. First, the K~ production cross section in near-threshold 
proton-proton reactions is expected [2] to be lower by over an order of mag- 
nitude than the K~^ cross section at the same distance from their respective 
thresholds. Second, antikaons {K" , K ) should have a higher absorption prob- 
ability in heavy ion collisions than kaons {K^ , K^)i primarily as a result of 
strangeness exchange reactions (e.g. K^N -^ Yn with F = A, S). Although 
charge exchange reactions (e.g. K'^n — * K^p) may occur, changes to the ob- 
served yield of a particular charge state will be largely compensated by the 
corresponding reverse process. In this context it has been shown in detailed 
model calculations that multi-step processes can not explain the enhanced K~ 
yield in heavy ion reactions [3]. A promising explanation of this discrepancy 
comes from various theoretical models [4] which suggest that, as a result of 
partial chiral symmetry restoration in a dense hadronic medium, antikaons are 
subject to strongly attractive, and kaons to slightly repulsive, forces. These 
effective interactions would lower the apparent K^ production threshold in a 
dense hadronic medium, and thus enhance the K~ yield in heavy ion collisions. 

However, the comparison of heavy ion with nucleon-nucleon collisions has 
relied on some assumptions because, in the latter case, the antikaon data 
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were either deduced from pp -^ ppK^K^ results, or taken as the sum of dif- 
ferent exclusive channels with one or two pions in the final state [2]. Fur- 
thermore, apart from a preliminary upper limit [5], there is a lack of pp an- 
tikaon production data in the regime of available energy most relevant in 
near-threshold heavy ion collisions (i.e. y^ — ^JsQ < 0.5 GeV ^ where ^/sq is 
the threshold of the particular channel under consideration) [6]. The result 
presented in this work represents the first cross section determination of the 
exclusive channel pp -^ ppK^K~ , and is equivalent to the inclusive K~ pro- 
duction cross section since no other purely hadronic channels including K~ 
are kinematically allowed at the beam momentum of this experiment (i.e. 
^ - ^ = 111 MeV < M^oc^). 



2 EXPERIMENTAL PROCEDURE 



2. 1 Apparatus 



A proton beam from the SATURNE proton synchrotron with momentum 
Vbeam = 3.67 GeV/c was directed onto a liquid hydrogen target of 2 cm 
length. All events with at least four charged particles in the final state were 
measured with the DISTO spectrometer, which is described in detail else- 
where [7]. Charged particles were tracked through a magnetic spectrometer 
which consisted of a dipole magnet (~ 1.0 T -m), two sets of scintillating fiber 
hodoscopes inside the field and 2 sets of multi-wire proportional chambers 
(MWPC) outside the field. Furthermore, two arrays of scintillator hodoscopes 
and water Cerenkov detectors were located behind the MWPCs. Particle iden- 
tification was performed using the correlation between the particle momenta 
and the Cerenkov light output (see also [7]). The large acceptance of the 
spectrometer (^ ±15° vertical, ~ ±48° horizontal in the laboratory refer- 
ence frame) guarantees a sizeable efficiency for coincident detection of four 
charged particles. The measurement of all particles in the final state allowed 
4-momentum conservation to be used in addition to particle identification for 
effective background suppression in order to identify the ppK^K~ final state. 



2.2 Data Selection 



Since the 4-momenta of all particles in the final state were measured, the 
events are kinematically over-determined. Therefore, 4-momentum conserva- 
tion can be used for a substantial background suppression by requiring that 
the proton-proton missing mass (M^^^^) be equal to the K^K~ invariant mass 
(M^f ). The distribution of (M^f )2 - {MZssY is plotted in Figure 1. The 



peak near (M^^)"^ — (M^^^)^ = results from events where the ppK~^K~ 
event hypothesis was correct, and is the basis for the cross section values 
quoted below. This peak is superimposed on a background resulting from im- 
perfect Ti — K separation in the Cerenkov detectors in a small fraction of 
events of the type pp — > pK^A -^ ppK^n^ or pp -^ ppn^ii^X. An estimate 
of the background underneath the ppK^K~ peak is given by the solid curve 
in Fig. 1. This estimate was determined by scaling the {M^^Y — (M^^^^)^ 
distribution, measured for all events before requiring kaon identification, by a 
factor 0.002 in order to match the data in Fig. 1 above 0.15 GeV'^/c^. From 
this histogram it is determined that the background accounts for about 13% 
of the yield in Fig. 1 with \{Mg^f - {MZssf\ < 0-09 GeV^/c^. This back- 
ground has been subtracted in the subsequent analysis as described below. To 
extract a total cross section from the data shown in Fig. 1, a correction for the 
detector acceptance and an absolute normalization must also be determined. 



2.3 Acceptance Corrections 



The correction of the measured yields for the detector acceptance has been 
evaluated by means of Monte Carlo simulations, which after digitization of 
the simulated detector hits, were processed through the same analysis chain 
as the measured data. The detector acceptance was determined as a multi- 
dimensional function of the relevant kinematic degrees of freedom of the parti- 
cles in the final state. After accounting for the azimuthal and reflection symme- 
tries, the detector acceptance was non-zero over the full kinematically allowed 
region. Thus, the acceptance correction is essentially independent of the actual 
phase space distribution of the final state, and has been determined assuming 
a uniform phase space distribution in the simulations. This has been verified 
by calculating the acceptance correction matrix using different initial distri- 
butions in the simulations, and observing that the measured yield varied less 
than the systematic error associated to the acceptance correction. Although 
eight linearly independent degrees of freedom are in principle required to fully 
describe the pp -^ ppX -^ ppK^K^ reaction, a five dimensional acceptance 
correction matrix is sufficient, in part because the cross section cannot depend 
upon the azimuthal orientation of the event. Furthermore, we have compared 
the observed angular distribution of the X -^ K^K~ decay with that from 
simulations. The comparison of data with simulations has a x^ /n = 1.4 for a 
S-wave distribution and a x^ /n = 17.6 for a P-wave distribution. Therefore, 
in the simulations uniform distributions with respect to these three angular 
variables were integrated over when determining the acceptance as a function 
of the remaining five variables. Finally, the raw data were corrected on an 
event-by-event basis, via a weighting factor determined from the simulated 
acceptance function for the appropriate kinematic bin. 



Since the acceptance varies as a function of the kinematic distribution of the 
final state, one can not simply reduce the acceptance to a single number with- 
out making assumptions about the phase space distribution of the particles 
in the final state. Nevertheless, a reasonable estimate of the average accep- 
tance can be determined for kaon pairs with an invariant mass equal to the 
phi meson mass (where most of the measured kaon pairs are observed) and an 
isotropic distribution in the other kinematic variables. In this case the product 
of the geometrical acceptance times the tracking reconstruction efficiency is 
21.4%. 



2.4 Background Subtraction 



The background contribution from non-ppK^ K~ events, shown as the solid 
curve in Fig. 1, must be subtracted from the data in order to determine 
the K^K~ yield. The subtraction was performed on the M^^ distribu- 
tion. For this, the M^^ distribution for the background events was deter- 
mined by applying the acceptance correction matrix (determined for simulated 
pp -^ ppK^K" reactions) to events subjected to the same kinematic require- 
ment \{M^^f - {M^i^,f\ < 0.09 GeV^/c^, but not to the kaon identification 
conditions with the Cerenkov detectors. The resulting distribution was scaled 
by the same factor (0.002) used by the solid curve in Fig. 1, and then sub- 
tracted from the acceptance corrected M^^ distribution that included both 
ppK^K" and background events. Finally, the M^^l^ spectrum, after full ac- 
ceptance corrections and background subtraction, is shown in Figure 2. The 
curves are fits to the data as described below. 



2.5 Absolute Normalization 



The absolute normalization of the ppK^K~ cross section was determined by 
measuring the yield relative to that of a simultaneously measured channel 
with known cross section. For this work the reference channel was the reac- 
tion pp — » pprj for which a large amount of data exist [8]. This method to 
determine the absolute normalization was chosen because it greatly reduced 
the large systematic uncertainty associated with the absolute calibrations of 
both beam intensity and absolute trigger efficiency. In order to provide the 
absolute cross section calibration the existing published data was first inter- 
polated to get the 7/ production cross-section at the beam momentum of the 
present measurement, and then the appropriate acceptance corrections were 
applied to determine the K^K^ /rj total cross section ratio from our yields. 

The cross section of the reaction pp -^ pprj has been interpolated to our beam 
momentum with several parameterizations that vary smoothly with beam mo- 



mentum. From these interpolations we estimate the exclusive pp -^ ppt] pro- 
duction cross section to be 135 ± 35/i6 at Pbeam = 3.67 GeV/c. The parame- 
terizations describe the existing rj production data well with the exception of 
the single measurement at pteam = 2.8 GeV/c by E. Pickup et al. [9] which is 
significantly underestimated. This discrepancy has been neglected since that 
measurement is subject to a large systematic error associated with a quite 
substantial background subtraction. 

In the present data, the f] meson has been identified in the M^^^ distri- 
bution for ppTT~^n~X events, after requiring that the four particle missing 
mass be consistent with M^^o as shown in Ref. [13]. The acceptance correc- 
tion of the pp -^ pprj -^ pp-K^-K^-K^ channel was performed similarly to the 
pp -^ ppK^K~ discussed above. The acceptance correction matrix for the pprj 
channel has four dimensions. These are sufficient to completely describe this 5 
body final state because the full set of 15 kinematic degrees of freedom (dof ) 
is reduced by four-momentum conservation (-4 dof), the azimuthal symmetry 
of the event (-1 dof), the requirement that M^^^ = M^ (-1 dof), the isotropic 
orientation of the rj meson decay plane (-3 dof), and the known matrix element 
(-2 dof) for the rj -^ 7r+7r^7r° decay [10]. 



3 Results 



3.1 Total K^K Cross Section 



After applying the full acceptance corrections to the pp -^ ppK^K~ and 
pp -^ pprj — > pp7i~^7r~7i^ channels, the K~^K~ /rj total cross section ratio is 
determined to be (1.5 ± 0.1{stat.) ± OA{sys.)) x 10^^. The systematic er- 
ror quoted here arises from the quadratic sum of the uncertainties on the 
rj and K^K^ background subtractions (15% and 5% respectively), relative 
acceptance correction (11%), trigger bias (10%), tracking efficiency (10%), 
and the Cerenkov particle identification efficiency (14%). This value has also 
been corrected for systematic bias effects due to the different scintillating fiber 
efficiencies for pions and kaons (-7.5%) and from the subtraction of events aris- 
ing from the target envelope (-1-5%). Based on the rj cross section estimated 
above in Sec. 2.5, a total cross section for the reaction pp -^ ppK^K~ of 
(0.20 ± 0.011 ± 0.08)/i6 has been determined, where the second error quoted 
here is the quadratic sum of the systematic uncertainty in the measured yield 
ratio and the absolute normalization uncertainty for the rj production cross 
section. 

The present total pp -^ ppK^K^ cross section value is plotted as the solid data 
point in Figure 3, where it is compared with estimates of K~ inclusive produc- 



tion cross sections at higher energy (open circles) taken from the hterature [2]. 
These additional data points have been deduced from other reactions either 
by assuming app^K-x = cr ^t^x' °^ ^^ taking the sum of exclusive channels 
with one or two pions in the final state. The solid curve, corresponding to a 
prediction from Sibirtsev et al. [11] using a one meson exchange model includ- 
ing vr, K, and K* mesons, accounts well for the K^ cross section measured 
in this experiment near the production threshold. In comparison, the total 
cross sections for inclusive K^ production shown as the open diamonds [12] 
in Fig. 3, are more than an order of magnitude larger at comparable distances 
above the K^ threshold. 



3.2 (f) Meson 



Near threshold the K~ meson is produced, to a large degree, by the decay of 
the meson as an intermediate state [13]. The fraction of the cross section 
from the resonant production can be determined from the M^^ distribution. 
The M^^ spectrum in Fig. 2 has been fit with the sum of a non-resonant 
contribution and a peak from the resonance. The shape of the non-resonant 
contribution was assumed to be given by the M^^ distribution for an ensemble 
of events that are uniformly distributed according to four body {ppK^ K') 
phase space. The shape of the resonance was given by the natural line- 
shape folded with a Gauss function to account for the detector resolution. 
When treating the width {a) of the Gaussian as a free parameter, we find an 
optimal fit with a = 3.5 ± 0.5 MeV/c^, in good agreement with simulations 
of the detector performance. 

The total K^K^ cross section value, as well as the resonant and non-resonant 
components, are summarized in Table 1. After correction for the corresponding 
partial width [14], the total meson production cross section is 0. 19 ± 0.014 ± 
0.08 /i6. The systematic errors quoted have the same meaning as explained 
above. 

Table 1 

Total exclusive production cross section for the reaction pp -^ ppK~^K~ at 3.67 

GeV/c and for the resonant (0 meson) and non-resonant contributions. 

Meson Species Cross Section [fib] 

total K+K- 0.20 ±0.011 ±0.08 

<j)^K+K- 0.09 ±0.007 ±0.04 

non-resonant K+K' 0.11 ± 0.009 ± 0.046 



4 SUMMARY 



In conclusion, the total cross section of the reaction pp -^ ppK^K~ has been 
determined for Pbeam = 3.67 GeV /c. This is the first cross section measure- 
ment of the pp -^ ppK^K~ channel near threshold, and is equivalent to 
the inclusive cross section at this beam momentum. The cross section de- 
termined here is more than a factor 20 lower than the measured [12] and 
calculated [15] pp —^ K^ + X cross section at the same CM energy above 
threshold. This large difference in the K^ and K^ meson production cross 
sections in proton-proton collisions is in strong contrast to the nearly equal 
cross sections measured in sub-threshold heavy ion collisions at the same dis- 
tance from the respective thresholds. Since this discrepancy between heavy 
ion collisions and proton-proton reactions has been interpreted as possible ev- 
idence for in-medium modifications of the K~-nucleon interaction, it would 
be very useful to study the evolution of the relative K^ and K~ meson yields 
at equal energies from the threshold versus increasing mass number of the 
colliding nuclei. 
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Fig. 1. Distribution of {M^^)^ — (M^^,^)^. The data points are for events with 
kaon identification and the sohd histogram is the scaled background deduced from 
events without kaon identification. 
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Fig. 2. Spectrum of M^^ after acceptance correction and background subtraction. 
The solid curve is a fit to the data with the sum of a non-resonant K'^K' contribu- 
tion (dashed curve) and the natural line-shape of the (j) resonance folded with the 
detector resolution. 
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Fig. 3. Total kaon and antikaon production cross sections as a function of the avail- 
able energy above the appropriate threshold. The present measurement is the solid 
data point and the open circles are antikaon points deduced from the literature. 
The curve is a model prediction described in the text and the open diamond points 
are cross section values for positive kaon production. 
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